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Abstract 

We present a renormalization group analysis to Einstein-Rosen waves or vacuum spacetimes 
with whole-cylinder symmetry. It is found that self-similar solutions appear as fixed points in the 
renormalization group transformation. These solutions correspond to the explosive gravitational 
waves and the collapsing gravitational waves at late times and early times, respectively. Based 
on the linear perturbation analysis of the self-similar solutions, we conclude that the self-similar 
evolution is stable as explosive gravitational waves under the condition of no incoming waves, 
while it is weakly unstable as collapsing gravitational waves. The result implies that self-similar 
solutions can describe the asymptotic behavior of more general solutions for exploding gravitational 
waves and thus extends the similarity hypothesis in general relativity from spherical symmetry to 
cylindrical symmetry. 
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I. INTRODUCTION 



Renormalization group analysis is a powerful tool for obtaining asymptotic behavior of 
solutions of partial differential equations jl|. Koike et al. j^, 3| have applied this method 
to a self-gravitating system in general relativity and successfully explained critical behavior 
together with the critical exponent in gravitational collapse. Ibanez and Jhingan 0, [sl 
have applied this method to inhomogeneous cosmology and analyzed the stability of scale- 



invariant asymptotic states (see also Refs. pi- 
Fixed points under renormalization group transformation are very important for the 
asymptotic analysis. The fixed-point solutions in general relativity generally correspond to 
self-similar solutions, which are defined as spacetimes that admit a homothetic Killin g ve ctor. 
Self-similar solutions arise from the scale-invariance of general relativity. Carr jol. llo| has 
conjectured that under certain physical circumstances, spherically symmetric solutions will 
naturally evolve to a self-similar form from complicated initial conditions, and this is termed 
as the similarity hypothesis. This conjecture has been strongly supported in the gravitational 
collapse of a perfect fluid by numerical simulation and linear stability analysis 11 1. 



Although the similarity hypothesis has been proposed originally in spherical symmetry, 
there is evidence that supports the validity of the conjecture also in cylindrical symmetry. 



Nakao et al. [12[ have numerically simulated the collapse of the dust cylinder and the sub- 
sequent emission of gravitational waves in special cylindrical symmetry. Their numerical 
simulation suggests that gravitational waves gradually approach a self-similar form. 

The vacuum Einstein equations in this symmetry can be analytically integrated, and 



solutions are called the Einstein- Rosen waves 



13| . Since this system admits two commutative 



spatial Killing vectors and still retains the dynamical degrees of freedom corresponding 
to gravitational waves, it provides us with a good toy model in which we can learn the 



physical nature of gravitational waves. Harada et al. have derived self-similar Einstein- 



Rosen waves and identified one of them with the attr actor solution that Nakao et al. [12 1 
reported. In this paper, we recover the self-similar Einstein-Rosen waves as a fixed point of 
the renormalization group analysis. Then, we perform the linearized perturbation analysis 
of these solutions (fixed points) and determine their stability analytically using an eigenvalue 
analysis. 

This paper is organized as follows. In Sec. II, we show the field equations for Einstein- 



Rosen waves. In Sec. Ill, we introduce the renormalization group analysis as a scaling 
transformation and self-similar Einstein-Rosen waves as fixed points. In Sec. IV, we apply 
linear perturbation theory to the self-similar solutions. In Sec. V, we analyze the behavior 
of the perturbations at the boundaries. In Sec. VI, we show that the self-similar Einstein- 
Rosen waves are stable and unstable at late times and at early times, respectively, under 
appropriate boundary conditions. In Sec. VII, we conclude the paper. We use the units in 
which c = 1. 



II. EINSTEIN-ROSEN WAVES 



We consider cylindrically symmetric spacetimes with the azimuthal Killing vector d/dcf) 
and the translational Killing vector d/dz. We additionally assume that these two Killing 
vectors are hypersurface orthogonal for whole-cylinder symmetry. Vacuum spacetimes with 



whole- cylinder symmetry are called Einstein-Rosen waves 
Rosen waves is given by 

ds"" = e^(^-^\-df + dx^) + e-^^xHct)^ + e^^dz\ 

The nontrivial Einstein equations take the following form: 

-ij,tt + i^.xx + -ij,x = 0, 

X 

l,x = X{lpl + 1pl), 

7, = 2xilj^^^lj,. 



13| . The line element in Einstein- 



(2.1) 

(2.2) 
(2.3) 
(2.4) 



We define new variables h{x,t) and f{x,t), which are more convenient to our analysis, by 



f = li-e2^ and h 



x^ 



respectively. Then, the line element is rewritten in the following form: 

ds^ = f{-df + dx^) + x{h^d(f)^ + h~^dz^), 

where / and h satisfy 



and 



f.t = ±2xf- 



hi = ±h 



h ... 



Ax'^ X f 



X f 



1/2 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



respectively, where (and hereafter), in the case of double sign, we should uniformly choose 
either an upper sign or lower sign. The wave equation (12 ■2p is identically satisfied from 
the above two equations. Hence, we have the desired evolution equations for the metric 
functions in a form suitable for a renormalization group analysis. 



III. RENORMALIZATION GROUP ANALYSIS 



Since we are interested in the solution in the asymptotic regime, we shall explore now 
the scale-invariant properties of the system. We consider the following scale transformation 
and define the scaled functions, /i(i)(x,t) and /(i)(x,t), as follows: 

X — )■ Lx, (3.1) 

t L% (3.2) 

h -> L^hiLx, L'^t) = /i(i)(x, t), (3.3) 

f^L'f{Lx,L^t) = f^L){x,t). (3.4) 

Here, L is the scaling parameter, and scaled quantities and /(l) satisfy the same 
evolution equations, i.e., Eqs. (12. 7p and (12. 8p . From the structure of dynamical system, it is 
easy to see that t scales in the same way as x, fixing a = 1. There is no further constraint 
on scaling exponents a and b. Putting t = 1 and redefining, in succession, Lx by x and L 
by t, in Eqs. O and ([33!), we find 

/i(x,t) = r"/i(L)(x/t, 1), (3.5) 
/(x,t)=t-V(L)(x/t,l) (3.6) 

in the scaling regime. Note that the quantities /i(i)(x/t, 1) and /(2,)(x/t, 1) are evaluated 
at some initial time {t = 1 here) and are determined as fixed points of the renormalization 
group equations as illustrated below. 

n 

In what follows, we adopt the mechanism developed by Bricmont et al. py]. Defining 
L = exp(r), we can rewrite the scaling behavior of Eqs. (13.30 and (13.41) into a set of coupled 
ordinary differential equations: 

dh( 



ah(^L) + xh(^L),x + h^L),t\^^^, (3.7) 



dr 

dr 
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bf(L) + xf(L),x + /(L),4=l ' (3-8) 



where the quantities on the right-hand side are evaluated at initial time and are functions 
only of X. The fixed-point solutions h*^^ and f"^^ are defined by 



dT 



dr 



(3.9) 



Therefore, we have fixed fixed points as solutions of a set of coupled nonlinear ordinary 
differential equations: 



a/l(L) + + = 0, 

&/(L) + a;/(L),. + 2x/(i)^vl/ = 0, 



h 



(3.10) 
(3.11) 



where \1/ is defined as 



^ = ± 



1 f{L),. 



X f(L) 



h 



(L),x 



h 



1 1/2 



(3.12) 



Here and henceforth, we have dropped superscript * from and f^i) for brevity. These 
equations decouple on recombination, and we have the following master equation for h(^Ly. 



h 



{L),x 



h 



X 



a^ + b-1/4: 



1/2 



(3.13) 



The other metric function /(i) can now be recovered as a solution to an ordinary differential 
equation: 

_ b 2{a' + b - 1/4) 2a^a^ + 6-1/4 
/(L) x(l-x2) + xVT^ ■ ^ ' 

We are interested in the domain < x < 1, which now corresponds to the spacetime 
inside the light cone < x < t. This restricts the parameter space to + 6 — 1/4 > 0. 
Equations f l3.13p and f l3.14p can be easily integrated to give the fixed points 



f[L) = CfX 



h^L)=x-^^'^{l-VT-x 



(3.15) 

(3.16) 



where c/ is a constant of integration and A = a"^ + b — 1/4. Recovering time dependence 
through Eqs. (13. 5 p and (13. 6p . we obtain 



h{x, t) = X-^^^it - Vt^ - X2)±^, 

f{x,t) = C^X-''+2A^±2aA(^2 _ ^2yA^^^ _ y^J3^)T2aA_ 



(3.17) 
(3.18) 



t ' 
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FIG. 1. The light cone spHts the spacetime into four parts. In the current paper, we focus on 
domains I and IV, i.e., the interior of the light cone. 

The light cone sphts the spacetime domain t > into two parts, and we have another solution 
corresponding to the region x > t > and A < 0. In fact. Fig. [T] shows that the light cone 
splits the spacetime into four parts if we consider both t > and t < 0. For the moment, 
we focus on domain I and later on IV, i.e., the interior of the light cone. For domains II and 
III, i.e., the exterior of the light cone, we have another solution corresponding to x > |t| > 
and A < 0, which is out of the scope of the present paper. 



These solutions coincide with the self-similar or homothetic solutions derived in Ref. jl^ . 
in which the solutions are further restricted by imposing the solution to be either with 
a regular axis or a conically singular axis. Then, the solutions are parametrized by two 
parameters k and A. The parameter k plays a more important role since it determines the 
physical nature of the spacetime, while A just controls the deficit angle of the axis. The 
correspondence between a, 6, and c/ in the above obtained expression and k and A in the 



expression given in Ref. IJ] is the following 



±A = K, a = K- 1/2, 6 = K, c) = 2-^^'e2^, (3.19) 

where we should note that A > 0. Note that we have started the analysis in the domain < 
X < t. Here, it should also be noted that both k = and 1/2 correspond to the flat spacetime, 
and K = —1/2 corresponds to the Kasner solution with exponents (—1/2,2/3,2/3) or the 
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locally rotationally symmetric Kasner solution. The x = t surface is a null singularity for 
< < 1/4 and 1/4 < < 3/8, a regular surface admitting an extension beyond it 
for 3/8 < < 1/2, and a null infinity for > 1/2. We can further divide the case of 
3/8 < < 1/2 into the case of = 1/2 - l/(4r2) {n = 2,3, ■ ■ ■) and otherwise. In the 
former, the x = t surface admits an analytic extension beyond it, while in the latter, this 
surface only admits a finitely differentiable extension. In general, these solutions describe 
exploding gravitational waves. If we flip the sign of t, these solutions describe the collapse 
of gravitational waves. The fully detailed description of the physical nature and causal 



structure of these solutions is given in Ref. 
of the numerical simulation by Nakao et al. 



14 1 . The parameter value which fits the result 
12| is K = —0.0206, for which the surface x = t 



corresponds to a null singularity, which is physically a shock of gravitational waves, and the 
spacetime cannot be extended beyond it. 

IV. LINEAR PERTURBATION ANALYSIS 

We introduce the perturbation of the fixed point solution as follows: 

\L) = hl){l + 5h^L)). (4.1) 

/(L) = /(L)(l + '5/{L)). (4.2) 

Assuming |5/i(L)| ^ 1 and ^ 1, we linearize the field equations for the perturbation. 

The linearized perturbed equations take the form 

/ 1 h*\ 
^h.+ [x + -^\5f,^, (4.4) 



dr V ^ h* ' 2x^ 



dSf 

2x 



1 /h*-' 



dr 

where we have dropped the subscript (L) from the metric functions for brevity. Note that 
all the quantities are evaluated at the fixed points (emphasized by an * ) and, hence, are 
known functions of x given by Eqs. ( 13. 15^ and (13. 16^ . 

We now compute the normal modes with the dependence 6h = e^'^ p{x) and 5f = e^'^a{x), 
where a; is a constant. From Eqs. (14.31) and (14. 4p . we obtain 

1 a \ 1 



Combining the two equations, we can solve for a as 

-2w(a T AVl ^ 2/\xVl 



(4.7) 



We can recover a master equation for p(x) from Eqs. (14.51) . (14. 6p and (14. 7p . i.e., 



P,xx ~^ 



1 u;(u; 

—P,x 



-p = 0. 



(4.8) 



x(l — x^) ' 1 — 

Again, the above is estimated at t = 1, and, hence, x should be regarded as x/t in the 
notation of Harada et al. [l^. Note that there is no a priori reason to assume that u is real. 
We define C = 1/ vT^-a?^ and g = C^p. Then, Eq. (14. 8 p transforms to 



d 



2^dg 



+ uj{uj + l)g = 0, 



(4.9) 



where < x < 1 corresponds to ^ > 1. This is the Legendre differential equation. Note 
that this equation has symmetry for the replacement of the parameter u O — (w + 1). 
The Legendre functions of the first and second kinds — -Pw(C) QcuiC): respectively — are 
solutions of the Legendre differential equation. General solutions can be expressed in terms 
of the Legendre functions simply as 



P 



2\uj/2 



x') 



ClPuj 



+ C2Qu 



1 



x^ 



(4.10) 



The above expression is convenient for Re(a;) > —1/2. For Re(ci;) < —1/2, due to the 
symmetry of the Legendre differential equation, we can express the solutions in terms of the 
Legendre functions with indices v = —{u + 1) in the range Re(z/) > —1/2. Then, we find 
the following expression that is more convenient: 



p=(l-a;2)-/2 



CiP_(^+l) 



+ C2Q^ 



1 



(4.11) 



The solution for another perturbation a of the metric functions can be constructed by 
Eq. dlTD. 



V. BEHAVIOR OF THE LINEAR PERTURBATIONS 

We impose the boundary conditions from the physical requirements and study the allowed 
range of u. To do this, we need the asymptotic behavior of Pu{z) and Qu{z), for Re(i^) > 
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—1/2 at z = 1 and z = oo, which are given in Refs. ISMlTj. The asymptotic behavior of 
Py(z) and Qu{z) at z = 1 are given by 



and 



Quiz) 



\ In ^ (z. ^ -1, -2, -3, . . . , ), 

d 



(5.1) 



(5.2) 



where 7 is the Euler number and V'(^) = 'dT^'^^^^^^ polygamma function. The 

asymptotic behavior of Pu{,z) and Qi,{z) for z — )■ cxd is given by 



P-1/2{Z) ^ 
P~l/2+ip{z) - 

where p 7^ 0, p G M and 



P^{z) ^ 2''ti~^I^T{v + 1/2)^71(1 + z/), (Re(z/) > -1/2) 
v^ln(8z) 



i{2z)~^l^ 



7rtanh(7rp) r(l + 2ip) 



r(iZ2+^p.)- - £(lZi^(2^)* 



r(l - 2ip) 



Q,{z) ~ 2-^^-^1/2^^^ ^ l)z-''-Vr(z/ + 3/2), 



(5.3) 
(5.4) 

(5.5) 
(5.6) 



respectively. 

To make the discussion clear, we hereafter concentrate on the stability of self-similar solu- 
tions with a regular or only conically singular axis. If we assume /t = for the parametriza- 
tion of Harada et al. [IJ], we find all the solutions of the family are fiat. Thus, we can assume 
K 7^ or A 7^ without the loss of generality. The asymptotic behavior of the solutions at 
X = and x = 1 are given as follows. 



A. Re(w) > -1/2 and cj / 



For Re(a;) > -1/2, we find 

p ~ ci + cshlnx + ln2 - 7 - V^(a; + 1)], (5.7) 
(J ~ -2ci(a^ A) + 2c2[-(aT A)[-lnx + ln2 - 7 - i){uj + 1)] ± A/cj] (5.8) 

at X = and 

r(a; + l) r(w + 3/2)^ ' 



a ^ -2cia2'^7r- 



T(u + 1) 



T{cu + 3/2) 



(5.9) 
.10) 
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at X = 1, where (and hereafter) only the relevant terms are shown. Note that the gamma 
function has no zeroes. 



B. Re(w) < -1/2 

For Re(w) < -1/2, we find 

p ^ Ci + C2[— Inx + ln2 — 7 — 'ip{—u)], (5.11) 
a ~ -2ci(a T A) + 2c2[-(a ^ A)[-lna; + ln2 - 7 - tpi-u)] ± A/u] (5.12) 

at X = and 

nr^c, 2~""^r(-a;-l/2) ^2nc.+i/2 , ^ 0^^+1/2 1/2 ^(-a;) 

. - -c,a2-^3/2^V2 r(-c.) , . 2— ^/M^ + l/2)r(-a;-l/2) _ 



C. cj = 

For = 0, from Eq. (14. 7p . we find p = const. Substituting this into Eq. (14. 6p . we find 
a = const. This corresponds to the rescaling of t and x by At and Ax, where A is a positive 
constant. Hence, this is not a physical mode but a gauge mode. In the following, we exclude 
this zero mode from the analysis for this reason. 



D. a; = -1/2 

For u = —1/2, we find 

p ~ ci + C2[-lnx + ln2 - 7 - ^(1/2)], (5.15) 

a ~ -2ci(a T A) + 2c2[-(a ^ A)[- Inx + ln2 - 7 - ^(1/2)] t 2A] (5.16) 

at a; = 0, and 

P- — Clin +—C2, 5.17) 

vr VI - V2 

2\/2 / 8 2A \ ^ , , 

cr~ ci a In — == =F ^== - V 27rc2a (5.18) 

TT V vl-a;2 Vl-a;^/ 

at a; = 1. 
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E. uj = 



-1/2 + ip {p^O and p G M) 



For w = — 1/2 + ip, we find 

p ~ Ci + C2[- Inx + ln2 - 7 - ^/^(l/2 + zp)], (5.19) 
a ~ -2ci(a T A) + 2c2[-(a ^ A)[- Inx + ln2 - 7 - ij{l/2 + zp)] ^ 2A] (5.20) 



at X = 0, and 



2-1/2+ip 



p ~ Ci«- 



r(i/2 + ^p)2^ _ r(i/2-^p)2 



7rtanh(7rp) [ r(l + 2?p) 



r(i - 2ip) 



r(i + ip) 



23/2— jp 



-1/2 + ip 



A 



ci r(i/2 + zp)2 ,^, r(i/2 + zp) 

+ r(i + .p) 



7rtanh(7rp) r(l + 22p) 



(5.21) 

:i - x')-^/2+i?^5.22) 



at a; = 1. 



VI. NORMAL MODES AND STABILITY ANALYSIS 

To see the stability in terms of normal modes, we impose boundary conditions at x = 
and X = 1 and see whether a growing mode exits or not. It is not trivial what boundary 
conditions we should impose on the solution. At least, the perturbation must be regular at 
both points; otherwise the linear perturbation scheme should break down. Noting a =]= A = 
— 1/2, the regularity at x = requires C2 = for Ke{uj) > —1/2 and C2 = for Re(cij) < —1/2 
because of the logarithmic divergence of the Legendre function of the second kind at x = 0. 
We additionally impose regularity at x = 1 and then find that only Ke{uj) > —1/2 is allowed. 
Conversely, for Re(a;) > — 1/2 and C2 = 0, the normal modes are regular and analytic both 
at X = and x = 1. 



A. Stability of the late-time solutions 



Our first motivation is to explain the numerical simulation by Nakao et al. 12|. The 
parameter value for the background is k = —0.0206, for which the surface x = t is an 
outgoing null singularity, as is shown by Harada et al. ij]. Since this null singularity is 
naked, we could inject an incoming wave mode there, in principle. However, the numerical 
simulation by Nakao et al. I2I will correspond to no injection of such an incoming wave. 



11 



This corresponds to the condition that p = o" = Oatx = l. This, as well as the regularity at 
X = 0, excludes all normal mode perturbations. This conclusion is valid also for all late-time 
(t > 0) solutions with any values of k under the condition of no incoming wave. Note also 
that if we do not impose the no incoming wave condition but only the regularity condition 
at X = 1, all normal mode perturbations with R,e{u)) > — 1/2 are allowed, and, hence, the 
self-similar evolution is strongly unstable. 



B. Stability of the early-time solutions 



So far, we have implicitly assumed late-time solutions, i.e., t > 0. However, it is also 
interesting to study the stability of early-time (t < 0) solutions in the context of gravitational 
collapse. The early-time solutions can be obtained by just flipping the sign of t. Since in this 
case X = —t is an ingoing null singularity or regular surface, we can impose the regularity 
at X = and x = 1. We have seen that under the regularity at both x = and x = 1, 
only Re(ci;) > —1/2 is allowed. Since the physical time evolution from t = — oo to t = 
corresponds to the decrease of r from r = oo to r = — oo, we can conclude that the self- 
similar early-time solutions are weakly unstable inside the null surface x = —t against the 
normal modes with eigenvalue u satisfying —1/2 < Re(u;) < 0. If Re(aj) > —1/2, all the 
values for u are allowed. However, during the time evolution, the most rapidly growing 
mode will dominate the perturbation. In this sense, growing modes with the growth rate 
Re(a;) = —1/2-1-0 will dominate the perturbation, and this suggests that the critical exponent 
is 2 if there appears a scaling law for the quantity of the mass dimension. However, there exist 
a countably infinite number of unstable modes; critical behavior in this system would be very 
different from those which have been studied so far. In the present context, we do not need 
to refer to the behavior of the perturbation outside the surface x = —t. This is consistent 
with the renormalizat ion group analysis for the spherically symmetric gravitational collapse, 
initiated by Koike et al. j2|, |3| . If we additionally impose the condition p = cr = Oatx = las 
in the late-time case, all normal mode perturbations are excluded, and, hence, the self-similar 
evolution becomes stable. 
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VII. CONCLUSION 



The numerical simulation by Nakao et al. [1^ strongly suggests that a self-similar grav- 
itational wave acts as an attractor in the vacuum region at late times after the explosive 
burst of cylindrically symmetric gravitational radiation in whole-cylinder symmetry. Moti- 
vated with this numerical result, we study the stability of self-similar Einstein- Rosen waves. 
There are exact solutions that describe the self-similar Einstein-Rosen waves. The param- 
eter values are identified for the numerical solution found in Ref. 13], in which the null 
surface x = t corresponds to a null singularity or a shock gravitational wave. 

We have analyzed the behavior of cylindrically symmetric linear perturbation around 
the self-similar solutions in terms of the normal mode analysis. We have found that the 
late-time self-similar solutions are stable inside the surface x = t oi gravitational waves at 
late times under the no incoming wave condition. This is the case for all the self-similar 
Einstein-Rosen waves at late times. 



We have also investigated the stability of self-similar Einstein- Rosen waves at early times, 
which describe the collapse of gravitational waves. This is important in the context of 
gravitational collapse. We find that self-similar Einstein-Rosen waves are weakly unstable 
inside the surface x = —t against regular cylindrically symmetric perturbations. 



The former case provides us with the demonstration that self-similar solutions that arise 
from the scale invariance of general relativity play an important role in the dynamics of 



gravitational waves and, hence, extends the applicability of the similarity hypothesis 10| . 



which was originally proposed for spherically symmetric spacetimes. We should note that the 
locally rotationally symmetric Kasner solution is known to be unstable against some sorts 
of homegenous but anisotropic perturbations [l8|. This suggests that the self-similarity 
hypothesis is a consequence of the restriction to the models with few enough degrees of 
freedom, which, hence, stabilize the self-similar solutions that are not stable in general. The 
renormalization group analysis together with linear perturbation scheme is quite useful to 
understand the asymptotic behavior of the somewhat general solutions of partial differential 
equations. 
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